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Abstract: This short article presents a mathematical formula required for metric 
corrections in image extraction and processing when using different length scale 

(N 

p | factors in three-dimensional space which is normally encountered in 

<L> 

C/} cryomicrotome image construction techniques. 
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In many scientific and industrial situations, the coordinates space is scaled by 

O 

different length factors in the three spatial directions, x, y and z; which affect the 
metric relations. For instance, in the cryomicrotomic image extraction techniques, 



t—{ the thickness of slices may be subject to errors or variations making the voxel size 

> 

in one direction larger or smaller than its standard size in the two other directions. 

o 

Consequently, the geometric parameters obtained from these images, which are 
based on the standard units of image space of an assumed cubic voxel unit, will be 
i— h contaminated with errors causing a distortion because of the missing scale factors 

• • 

• 1-h required by the isotropy of the physical space. 

>< 

In the following we present a simple case based on real-life cryomicrotomic 
image construction algorithms in biomedical applications where vasculature trees 
are obtained by computing the radius of each vessel in a number of rotational 
steps through a whole circle and the results are then averaged to obtain the final 
radius [1, 2]. As these rotational steps are oriented differently in the 3D space, the 
contribution of the length scale factors will vary from one orientation to the other 
and hence a scaling correction is required to obtain the correct radius. 
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There are several possible ways for deriving a formula for this correction; these 
include the use of rotational matrices, and circle projection on the three standard 
planes (i.e. xy, yz and zx) to obtain three ellipses from which the three spatial 
components can be computed. However, an easier and more efficient way is to find 
a parameterized form of an intersection circle between a plane perpendicular to 
the vessel axis and the vessel itself, which in essence is equivalent to a great circle 
intersection of this plane and a sphere having the same radius as the vessel [3]. 
This method of derivation is outlined below. 

For a regular cylindrical straight vessel oriented arbitrarily in 3D space and 
defined by its two end points P±(xi,yi, z±) and ^2(^2, 2/2, £2) and radius r, a free 
vector oriented in its axial direction is given by 



while a plane perpendicular to this vector and (for simplicity with no loss of 
generality) passing through the origin is given by 



Now, in 3D space a parameterized circle of radius r centered (with no loss of 
generality) at the origin and lying in a plane identified by two orthonormal vectors 
b and c is given by the equation 




(1) 



a x x + a y y + a z z = 



(2) 



r [cos(t)b + sin(i)c] 



< t < 2tt 



(3) 



that is 
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(r [cos(t)b x + sm(t)c x ] , r [cos(t)b y + sin(t)c ?/ ] , r [cos(t)6 2 + sin(t)c z ]) < t < 2n 

(4) 

To find b and c, a formal orthogonalization process, such as Gram-Schmidt, 
with normalization can be followed where random vectors non-collinear to a can 
be used. However a more convenient way is to find an arbitrary non-trivial vector 
lying in the plane by inserting arbitrary values for two variables (e.g. x = 1 and 
y — 1 ) in the plane equation and solving for the other variable (z) followed by 
normalizing through the division by its norm. If this vector is considered b, then 
vector c is found by taking the cross product a x b and normalizing. 

If the following length scale factors: a, (3 and 7 are introduced on the x, y and z 
directions respectively, then the distorted radius, r', at a random orientation t — 6 
is given by 

r' = ryj (a [cos(9)b x + sm(0)c x }) 2 + (0 [cos^)^ + sm(0)c y }) 2 + (7 [cos(9)b z + sm(9)c z }) 2 

(5) 

and hence the actual radius, r, is given by 



yj (a [cos(9)b x + sin(fl)cj) 2 + {(3 [cos(9)b y + sm(9)c y }) 2 + (7 [cos(#)6 2 + sin(#)c z ]) 2 

(6) 

As the image construction algorithm computes r' at iV rotational steps (e.g. 
360 steps corresponding to 360°) and averages the results, to restore the corrected 
radius r, this correction should be introduced at each one of these steps. In an 
ideal situation where rotational symmetry holds, only one quarter of these steps, 
^, is required, resulting in a substantial computational economy. However due to 
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the measurements and algorithmic errors at each step, it may be safer to maintain 
the N steps as the errors are expected to level out or diminish by applying this 
process through the whole circle. 

This correction can also be extended to use for post processing correction by 
applying the correction on the final averaged radius following a correction-free 
extraction process. For N rotational steps we have 

N N 

r 'i = r Yl V ( Q Mfc) 6 * + sin(^)cj) 2 + (/3 [cos^i)^ + sin(^) Cj/ ]) 2 + (7 [cos(^)^ + sin(^)c 2 ]) 2 

i i 

(7) 

Since the averaged post processing radius is 

Ef K 

Rav - ^ (8) 

the actual radius is then given by 



TV / 2 2 2 

J2i y ( a [cos(^)& x + sin(^)c x ]) + (f3 [cos(9i)b y + sm(9i)c y }) + (7 [cos(6»i)6 2 + sin(^)c z ]) 

(9) 

Although post processing correction may not result in computational efficiency 
it may be more convenient and useful to use when the non-corrected data are 
already obtained with no requirement to repeat the extraction process. 

It should be remarked that this correction can be applied in general to correct 
for this type of distortion regardless of the number of steps (single or multiple) 
and the shape of the object as long as the x, y and z components of the position 
vector can be obtained for each point in space required to trace the path of the 
distorted shape. This process can also be extended from discrete to continuous by 
substituting the summations with integrations with some other minor modifications 
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to account for this correction in analytical contexts rather than numerical discrete 
processes. 



6 



References 

[1] R.D. ter Wee, H. Schulten, M.J. Post, J.A.E. Spaan. Localization and 
visualization of collateral vessels by means of an imaging cryomicrotome. Vascular 
Pharmacology, 45(3): e63-e64, 2006. 

[2] B. Bracegirdle. A History of Microtechnique: The Evolution of the Mi- 
crotome and the Development of Tissue Preparation. Science Heritage Ltd, 2nd 
edition, 1986. 

[3] G.B. Thomas, R.L. Finney. Calculus and Analytic Geometry, Addison Wes- 
ley, 9th Edition, 1995. 



